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Построены математическая и имитационная модели циркадианного ритма систолического артериального 
давления, а также стохастическая и имитационная модели регуляции систолического артериального давле- 
ния в моменты стрессовых ситуаций. 
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Введение. Гипертоническая болезнь (ГБ) является актуальной проблемой кардиологии. В России 
около 40 млн человек страдают этим заболеванием. Одной из основных причин смертности тру- 
доспособного населения страны являются осложнения ГБ (инфаркт миокарда, хроническая сер- 
дечная и почечная недостаточность, атеросклероз и др.). В связи с тем, что существенное изме- 
нение артериального давления (АД) способствует развитию таких осложнений, актуальным явля- 
ется построение математических и компьютерных имитационных моделей регуляции АД у лиц, 
страдающих ГБ. 

В большинстве работ математическое описание медико-биологических процессов дано в 
терминах обыкновенных дифференциальных уравнений [1]. В этом случае рассматривается пове- 
дение объекта без учета случайных воздействий внешней среды. Организм человека - это дина- 
мическая система с элементами хаотичной структуры, в которой невозможно со стопроцентной 
вероятностью предсказать ее состояние в каждый момент времени. Поэтому детерминированные 
математические модели не могут адекватно описать ее поведение. При изучении динамики АД 
особый интерес представляет имитационное моделирование на основе стохастических диффе- 
ренциальных уравнений [2]. 

Цель исследования — построение стохастической и компьютерной имитационной моделей 

регуляции систолического АД (САД) в моменты стрессовых ситуаций на основе эксперименталь- 
ных данных. 
Объекты исследования. В исследовании принимали участие 71 женщина и 102 мужчины. По 
результатам суточного мониторирования АД (СМАД) пациенты были разделены на 6 групп без 
учета возраста: лица без ГБ (25 женщин и 31 мужчина), лица с ГБ без терапии (32 женщины 
и 56 мужчин) и лица с ГБ на терапии (14 женщин и 15 мужчин). 

СМАД проводили с использованием монитора «ВРЁаб МнСДП-3» (ООО «Петр Телегин», 
Н. Новгород). Пациенты наблюдались в Областном центре артериальной гипертонии Ульяновско- 
го клинического госпиталя ветеранов войн в 2008-2010 гг. 

Для статистической обработки данных использовали пакеты прикладных программ ВРЁаБ 

у. 3.0 (ООО «Петр Телегин», Н. Новгород), 5ТАТТЬТТСА 8.0 (Зао(, Тпс., Ч5А). Для реализации 
компьютерной модели был выбран язык программирования Войапа Веры 7.0 (Воцапа боЁмаге 
СогрогаНоп, Ц5А). 
Статистический анализ показателей СМАД. У лиц без заболеваний сердечно-сосудистой сис- 
темы суточный профиль АД характеризуется циркадианными колебаниями с периодом (24-4) ча- 
са и двухфазным ритмом «день-ночь». При этом степень ночного снижения САД не менее 10 %. 
Проведенные исследования показали, что суточный ритм АД существенно изменяется при ГБ. 





* Работа выполнена в рамках федеральной целевой программы «Научные и научно-педагогические кадры инновационной 
России» на 2009-2013 гг., а также при поддержке Министерства образования и науки РФ в рамках постановления прави- 
тельства РФ № 218. 
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В норме у лиц без ГБ и лиц с ГБ на терапии гомеостатические механизмы должны поддер- 
живать значения АД в определенных границах в течение некоторого периода, который мы назо- 
вем промежутком стабилизации АД. У таких пациентов предполагается наличие дневного и ноч- 
ного промежутков стабилизации САД и диастолического АД (ДАД). 

У лиц с ГБ без терапии механизмы гомеостаза АД нарушены, поэтому давление значи- 
тельно колеблется, что приводит к структурным изменениям его суточного профиля и отсутствию 
промежутков стабилизации или уменьшению их длины. 

Были проанализированы суточные профили САД и ДАД на наличие периодов стабилиза- 
ции. Большинство пациентов, у которых есть хотя бы один промежуток стабилизации САД, 
оказалось в группе женщин с ГБ на терапии (42,9 %) и группе лиц без ГБ (36,0 % женщин и 
41,9 % мужчин). Меньше всего таких пациентов было в группе женщин с ГБ без терапии (18,7 %) 
и группе мужчин с ГБ на терапии (20,0 %). 

Пациентов, у которых есть промежутки стабилизации ДАД, больше всего было в группах 
женщин без ГБ (60,0 %), женщин и мужчин с ГБ на терапии (71,4 % и 60,0 % соответственно). 
Меньше всего таких пациентов было в группе мужчин без ГБ (38,7 %). Согласно полученным ре- 
зультатам, почти в каждом случае у ДАД больше периодов стабилизации, чем у САД. 

Нами проанализирована статистическая взаимосвязь между длинами промежутков стаби- 
лизации АД и другими гемодинамическими показателями методом ранговой корреляции по Спир- 
мену (критический уровень значимости р < 0,24). Согласно исследованию, чаще всего между со- 
бой коррелировали: длины дневных промежутков стабилизации САД и ДАД; длина промежутка 
стабилизации ДАД и вариабельность ДАД (за дневной и ночной периоды суток соответственно). 

Вариабельность АД представляет собой случайные колебания, которые невозможно пред- 
сказать на 100 %. У пациентов с ГБ высокая вариабельность АД увеличивает риск сердечно- 
сосудистых заболеваний. В каждой группе был проведен статистический анализ значений вариа- 
бельности АД. Больше всего пациентов с нормальной вариабельностью было у лиц без ГБ (80,0 % 
женщин и 90,3 % мужчин). 


Математическая модель циркадианного ритма САД. Пусть случайный процесс У = (у ) 


1 


описывает динамику САД в мм рут. ст., время # измеряется в часах. Выделим у процесса У детер- 
минированную и стохастическую составляющие. Циркадианный ритм САД опишем детерминиро- 
ванным процессом С(1). Тогда случайные изменения САД (вариабельность) определим с помощью 
стохастического процесса. 

На практике для анализа колебаний САД традиционно используется математический ме- 
тод, который заключается в аппроксимации циркадианного ритма САД синусоидой [3]. Однако 
при обработке экспериментальных данных можно заметить, что график околосуточного ритма 
САД значительно отличается от графика синусоиды. В нашем исследовании рассмотрен новый 
подход к математическому моделированию динамики САД с учетом особенностей его суточного 
профиля. 

В качестве приближения циркадианного ритма САД воспользуемся непрерывной функцией 
С(=ЛАз(И)+р с периодом Т, где Т - время, в течение которого проводилось СМАД, ^.> 0 — 
амплитуда колебаний САД, р > 0 -— средний уровень САД по данным суточного мониторирования 
на отрезке [0; 7]. 

Пусть даны две числовые последовательности {11 и {у}, г меняется от 0 до п — 1. Число 
у; - значение САД, полученное в момент времени 1, е [0; 7] для любого г. Тогда параметры моде- 





У.=Р|. 


1 п-1 
ли: = 2 , ^ = тах, 
Непрерывная функция э(1) имеет вид: 
а, зт($()),0<1<1°, 
Е Л 
а, зш(()}1 т. 
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где промежутки [0; #] и (2; Т| рассматриваются как дневной и ночной периоды суток соответст- 
венно, {е (0; Т). Параметры а», { = 1, 2, определяются по методу наименьших квадратов (МНК). 
Для типов ночного снижения САД «поп-Фррег», «Чррег» и «оуег-@ррег» а; > 0, 1=1, 2. Для типа 
ночного снижения САД «паН-реаКег» а! > 0, а> < 0. Непрерывная неубывающая функция $Ф(®), 
агсзш(мм рт. ст.) — фаза САД, Ф(® е [0; 2] при ге [0; Т|. 

Определим ф(р) следующим образом. Если по данным мониторирования на графике суточ- 
ного профиля САД: 
пк 


1) нет промежутка стабилизации при 1 е [0; #], то Ф(#)=АЕ, где параметр А, =-, 


= 


агсзт(мм рт. ст.)/ч; 
2) есть промежуток стабилизации при ге [0; #], то 


| 
бы 


(^-Ф,)1 фт! 2 0 

? ее 

Здесь [1; 2] - дневной промежуток стабилизации САД. Параметр с! найден как выбороч- 

УР 
^, 





ное среднее значение последовательности {5;}, где 5, = для любого 1. При этом все 5; рас- 


: С. С. 
сматриваются в моменты времени & е [1; Р]. Параметр Ф. сии В случае, если —>1, 
а, а 


примем а! = 1; 
3) нет промежутка стабилизации при { е (#; Т|, то (+) =^1-+Ь,, где параметры 
п(Т-2г) 
Ве" 
4) есть промежуток стабилизации при 2 е и Т|, то 
(Ф;-п)! ле-фг 
3 0 + 3 0 
Е Е 
Ф(Р)=4Фь,Г <<, 
(2^-9,)'  ФТ-2ли 
та ТЫ 
Здесь [2; #] - ночной промежуток стабилизации САД. Параметр с› найден как выборочное 
среднее значение последовательности {5;}, при этом все 5; рассматриваются в моменты времени 


К Я агсз(мм рт. ст.)/ч; Ь, = агсзш(мм рт. ст.); 





ева. 





И, 


3. А 5 | ©2 С в 
1: Е [Г; []. Параметр ф, =п- агсзш| — |. Если — <-1, примем а, =1 для типов ночного сни- 
а, а> 


жения САД «поп-@ррег», «@ррег» и «оуег-@ррег» и а, =-—1 для типа ночного снижения САД 


«па -реакКег». 
Функция Ф(0 е [0; л] приге [0; #], Ф(0) Е (п; 2] приге (#; Т| и является непрерывной на 
отрезке [0; 7]. 
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При построении циркадианного ритма САД С(Ф) на любом промежутке времени [0; К], 


ат (Фа (#)),ЛТ << +, 


К > Т, функция &,(!)= фаза САД ф,(г)=$(1)+2л/ при 


а, т (Фа (1), г + /Т <1<(1+0Т, 
ге ПГ; О-ЮТ, где 1=0,т, "|= (функция [1] - целая часть числа Г). Тогда 


С() =, (К +р. 

Математическая модель регуляции САД. Часто причиной повышенной вариабельности АД 
являются стрессовые воздействия на организм. Стресс приводит к усиленной секреции в мозговом 
веществе надпочечников «стрессового» гормона адреналина, который затем выбрасывается в 
кровь, повышая САД. 


Пусть случайный процесс Х =(Х,)_, — концентрация адреналина в крови в момент вре- 


1>0 
мени + нмоль/л. Процесс выброса надпочечниками адреналина в момент стресса имеет вид пал, 
где л= (л,)_,- стандартный пуассоновский процесс с интенсивностью и > 0, п > 0 — коэффициент 


роста, нмоль/л. Тогда (—с)Х', — гомеостатические механизмы, понижающие концентрацию адрена- 
лина в крови после стрессовой ситуации, @ > 0 — коэффициент отрицательной обратной связи 
(коэффициент затухания). 

Процесс повышения САД под воздействием адреналина имеет вид ВХ,, В > 0 — коэффици- 
ент роста, (мм рт. ст.)-л/нмоль. Чтобы нормализовать давление, организм запускает гомеостати- 
ческий процесс —(ЙУ, который частично компенсирует стрессовое воздействие, где 

у» /Т << +, В | . ы 
(2) = т: 1=0,т, у;> 0, #=1, 2, = коэффициенты отрицательной обратной 

у ВЕЗЕТ 
связи (коэффициенты затухания), регулирующие скорость, с которой нормализуется давление в 
течение дневного (1= 1) и ночного (1 = 2) периодов суток. 

Другие факторы внутренней и внешней среды, влияющие на уровень САД (биологически 
активные вещества и гормоны, физические нагрузки, качество питания, метеоусловия, гиподина- 


9,/Г<#<Ё +, 
мия ит. д.), обозначим как о<(ВУДИ’, где с(!)= : 1=0,т, о; 0, #=1,2, = 
ОЕ РЕ ЗЕОИТ, 
пропорциональные коэффициенты, регулирующие изменение вариабельности САД при переходе 
от дневного (7 = 1) к ночному (7 = 2) периоду суток, а И =(/,)_, — стандартный винеровский 





1>0 
процесс. Процессы л и И’ независимы. 
Математическая модель регуляции САД в моменты стрессовых ситуаций представляет со- 
бой систему стохастических дифференциальных уравнений 
АХ, = пал, — @Х ‚4, 
ау =ВХ.@-у(1)У4+ с (г) Уа +С( 


с начальными условиями Ху > 0, > 0. 

Предполагается, что в данной системе процесс Х ненаблюдаемый, а наблюдать можно 
только процесс У, который содержит неполную информацию о процессе Х. Такая система называ- 
ется частично наблюдаемой. 

В предложенной модели случайный процесс Х не является непрерывным. Он совершает 
скачки в моменты времени, совпадающие с моментами скачков пуассоновского процесса л. Про- 
цесс У непрерывен, имеет множественные разладки. Разладка наблюдаемого случайного процесса — 
это изменение его характеристик (математического ожидания, дисперсии). Нами под разладкой 
понимается воздействие стресса на организм. Тогда моменты возникновения разладок - это мо- 
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менты стрессовых ситуаций. У процесса У они совпадают с моментами скачков пуассоновского 
процесса х. 

На рисунке изображен график траектории процесса У. В данном примере показатели 
СМАД взяты у пациента из группы женщин с ГБ на терапии, тип ночного снижения САД «аррег». 
Имеются два промежутка стабилизации САД. Общие параметры модели: К = Т = 19,4 ч, дискрет- 
ность по времени Л, = 0,01 ч. Параметры циркадианного ритма САД С(®: ^, = 37,16 (безразмерная 
величина), р = 156,16 мм рт. ст.; Ё = 9,8 ч, 1 = 0,67 ч, Р= 8,72 ч, Ё = 11,38 ч, Ё = 15,38 ч; 
К = 0,32 агсят (мм рт. ст.)/ч, а, = 0,49 (безразмерная величина), с, = 0,36 мм рт. ст.; 
А = 0,33 агсят (мм рт. ст.)/ч, 5> = —0,055 агсзш (мм рт. ст.), а> = 0,75 (безразмерная величина), 
с› = -0,65 мм рт. ст. Параметры процесса пл: и = 3 (безразмерная величина), ц = 2 нмоль/л. 
Параметры процесса Х: Х = 0, а = 10 (безразмерная величина). Параметры процесса У: 
У = 157 мм рт. ст., В = 10 (мм рт. ст.)-л/нмоль, безразмерные величины у, = 1,02, \› = 1,04, 
с! = 0,035, о> = 0,075. Коэффициенты цв, т, а, В, у1, 2, ол, > были найдены с помощью МНК. 


САД, ммрт. ст. 1 
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Траектория процесса У 


Проанализируем влияние параметров циркадианного ритма САД С(р на динамику траек- 
торий процесса У и его математического ожидания. Параметры модели А, А› являются вспомога- 
тельными. Они необходимы для того, чтобы функция =(1) проходила через заданные точки (0; 0), 
(6; 0), (Т. 0). Параметры аи, а — амплитуды колебаний функции (0) в течение дневного и ночного 
периодов суток, т.е. на промежутках [0; 1] и (6; Т| соответственно. Параметр /, — пропорциональ- 
ный коэффициент, определяющий амплитуду колебаний циркадианного ритма САД С(®) на всем 
промежутке [0; 7]. Чем больше абсолютные значения аи, а> и ^,, тем сильнее колебания циркади- 
анного ритма С(0) и тем больше отклонение математического ожидания процесса У от прямой 
Ки) = р. При построении компьютерной имитационной модели математическое ожидание процесса У 
было определено как функция, являющаяся средним арифметическим большого числа траекторий 
(например, одной тысячи) данного процесса. Параметры ст, с› необходимы для того, чтобы задать 
на графике траектории процесса У дневной и ночной промежутков стабилизации САД в течение 
периодов времени [1'; Р] и [Р; ] соответственно. Параметр р — выборочное среднее эксперимен- 
тальных значений САД на отрезке [0; 7]. Чем больше величина р, тем большие значения прини- 
мает математическое ожидание процесса У. 

При увеличении ц, интенсивности пуассоновского процесса х, увеличивается число скач- 
ков у процессов л и Х. Это приводит к увеличению числа разладок на графиках траекторий про- 
цесса У, при этом увеличиваются значения процесса У. Параметры п, В — коэффициенты положи- 
тельной обратной связи, при их увеличении значения процесса У возрастают. Параметры а. у, у› — 
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коэффициенты отрицательной обратной связи, при их увеличении значения процесса У умень- 
шаются. Параметры о1, с› характеризуют случайные воздействия внутренней и внешней среды на 
динамику процесса У. Чем больше абсолютные значения ст, с, тем сильнее траектории процесса 
У отклоняются от его математического ожидания. 

Выводы о характере влияния параметров модели на поведение траектории процесса У 
нашли подтверждение в ходе компьютерного имитационного моделирования. 

Предложенная нами математическая модель отличается от большинства моделей, описы- 

вающих процесс регуляции САД. Она учитывает такие особенности суточного профиля САД, как 
промежутки стабилизации, степень ночного снижения, изменение вариабельности в течение су- 
ток. Динамика САД задается с помощью непрерывного случайного процесса, имеющего множест- 
венные разладки с частичной компенсацией, в то время как процессы с разладками редко исполь- 
зуются при описании медико-биологических объектов. 
Заключение. Практическая значимость исследований состоит в построении математической и 
компьютерной имитационной моделей регуляции САД в моменты стрессового воздействия на ор- 
ганизм. Модели могут быть усложнены в зависимости от конкретной решаемой задачи, например, 
при добавлении в систему новых объектов, формирующих суточный профиль САД. 
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